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We pursue applications of the light-front reduction of current matrix elements in the Bethe- 
Salpeter formalism. The normalization of the reduced wave function is derived from the covariant 
framework and related to non- valence probabilities using familiar Fock space projection operators. 
Using a simple model, we obtain expressions for generalized parton distributions that are continuous. 
The non-vanishing of these distributions at the crossover between kinematic regimes (where the 
plus component of the struck quark's momentum is equal to the plus component of the momentum 
transfer) is tied to higher Fock components. Moreover continuity holds due to relations between Fock 
components at vanishing plus momentum. Lastly we apply the light-front reduction to time-like 
form factors and derive expressions for the generalized distribution amplitudes in this model. 
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I. INTRODUCTION 

More than a half century ago, Dirac's paper on the 
forms of relativistic dynamics jlj] introduced the front- 
form Hamiltonian approach. Applications to quantum 
mechanics and field theory were overlooked at the time 
due to the appearance of covariant perturbation theory. 
The reemergence of front-form dynamics was largely mo- 
tivated by simplicity as well as physicality. The light- 
front approach has the largest stability group j| of any 
Hamiltonian theory. Today the physical connection to 
light- front dynamics is transparent: hard scattering pro- 
cesses probe a light-cone correlation of the fields. Not 
surprisingly, then, many perturbative QCD applications 
can be treated on the light front, see e.g. ||. Outside this 
realm, physics on the light cone has been extensively de- 
veloped for non-perturbative QCD Gj as well as applied 
to nuclear physics ||, E|. 

Recently we investigated current matrix elements in 
the light-front Bethe-Salpeter formalism |7j , at task orig- 
inally attempted in As an example, we considered 
the ladder approximation for the covariant kernel in the 
weak-binding limit. We calculated (1 + l)-dimensional 
electromagnetic form factors in this reduction scheme to 
demonstrate the replacement of non-wave function ver- 
tices (as coined in Q) with contributions from higher 
Fock states. These higher Fock states originate from the 
light-front energy pole-structure of the Bethe-Salpeter 
and photon vertices. 

Below, we take the model into (3 + 1) dimensions and 
make clear the connection to higher Fock components by 
explicitly constructing the three-body wave function from 
our expressions for form factors. Additionally we show 
how the normalization of the covariant Bethe-Salpeter 
equation turns into a familiar many-body normalization 
in the light-front reduction. The immediate application 
of our development for form factors (in a frame where the 
plus-component of the momentum transfer is non-zero) 
is to compute generalized parton distributions. Connec- 
tion is again made to the Fock space representation and 



continuity of the distributions is put under scrutiny. The 
formalism is also employed to obtain time-like form fac- 
tors. This latter application in interesting since no Fock 
space expansion in terms of bound states is possible. 

The paper is organized as follows. First in section || we 
review the normalization of the covariant Bethe-Salpeter 
wave function and then proceed to derive the light-front 
reduced version. This necessitates a review of the light- 
front reduction notation introduced in [To[| . The normal- 
ization condition resembles a diagonal matrix element of 
a pseudo current and has contributions from higher Fock 
states. We calculate the explicit normalization condition 
for the ladder model at next-to-leading order in pertur- 
bation theory. The connection to the familiar many-body 
Fock state normalization is made in Appendix Next 
in section III, we use the reduction scheme to calculate 
generalized parton distributions. We do so by using the 
integrand of the electromagnetic form factor calculated 
in an arbitrary frame (these expressions in (3 + 1) dimen- 
sions are collected in Appendix A] where the leading-order 
bound-state equation for the wave function also appears) . 
We discuss the continuity of these distributions in terms 
of relations between Fock components at vanishing plus 
momentum. Connection is also made to the overlap rep- 
resentation of generalized parton distributions. Appli- 
cation of the reduction scheme to time-like form factors 
and their related generalized distribution amplitudes is 
presented in section IVl Finally we conclude with a brief 
summary (section |V ) . 



II. NORMALIZATION 

In [ |l~0| | the relation between the four-dimensional 
Bethe-Salpeter wave function and the reduced light-front 
wave function was presented. The normalization of the 
reduced wave function, however, was not discussed in 
much detail and will be addressed below. 

Let G{R) denote the two-particle disconnected prop- 
agator, where R labels the total momentum. Thus be- 
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tween effective singie-particle states of momenta p and 
k, we have (p\G(R)\k) = (2tt) 4 S 4 { p - k)G(k, R), where 
G(k, R) = d(k)d(R — k) and the scalar single-particle 
propagator is 



d(k) 



(k 2 - m 2 )[l + (k 2 - m 2 )f(k 2 )} + 



if 



(1) 



Here the function f(k 2 ) characterizes the renormalized, 
one-particle irreducible self-interactions and for simplic- 
ity shall be ignored below. The Bethe-Salpeter equation 
for the bound-state amplitude \*£ R ) with mass R 2 = M 2 
reads 



= G(R)V(R)\* R ) 



(2) 



where V(R) is the irreducible two-to-two scattering ker- 
nel (which we shall often call the potential). 

From the behavior of the reducible four-point function 
near the bound-state pole R 2 = M 2 one can deduce the 
covariant normalization condition by application of 
l'Hopital's rule 



(3) 



R 2 =M 2 



The normalization (|^) takes the form of a diagonal matrix 
element of a pseudo current. 

The light-front reduction is performed by integrat- 
ing out the minus-momentum 1 dependence with the 
help of an auxiliary Green's function G(R). For sim- 
plicity, we denote the integration J 4^(k~ \0(R) = 

O(R). With this notation, we will always work in 
(3 + l)-dimcnsional momentum space for which the 
only sensible matrix elements of O(R) are of the form 

(k+,k- L | 0(R)\p-,p + ,p ± ). The operator O(R) is de- 
fined similarly. To obtain light-front time-ordered per- 
turbation theory one chooses 



G{R) = G(R) g~\R) G(R), 



(4) 



where g(R) = G(R) . This form of G allows for a sys- 
tematic approximation scheme for the light-front energy 
poles of the Bethe-Salpeter vertex . Lastly one defines 
an auxiliary kernel W(R) by jl2| 



W(R) = V(R) + (G(R) - G[R))w(R). 



(5) 



In what follows we shall omit total four-momentum 
labels since they are all identically R. The normalization 



For any vector a M , we define the light-cone variables = 
^(a°± a 3). 



condition for the reduced wave function is then deduced 
by using the conversion |Tc|| 



(6) 



1+ [G— G )W \G \ 1R ) 



and the definition of the reduced wave function \ipR), 
namely 



\* R )=g(R)\ lR ). 



(7) 



Hence taking the plus component of Eq. 

2tR + = (7r\ \g^1 + W(G-G)\ 
( d 



dR- 



G^-V 



1 + (G-G)WJG | 7 .R>. (8) 



The complicated normalization condition is indicative of 
the effects of higher Fock space components. To see this 
explicitly, we work in the ladder model in perturbation 
theory for which 



V(k,p) 



(k — p) 2 — /i 2 + ie ' 



(9) 



Notice dV/dR^ = 0. 

Let us start with the contribution at leading order in 
G — G to the reduced wave function's normalization. 



;<7k| \g( 



d 



■G- 1 ]G 



2R+" 1 \dR 
To perform the integration, we note 

d 



\1R) = 1. 



(10) 



-G- 1 {k,R) = -2iR + d- 1 (k){l-x) 1 (11) 



dR 

where we have customarily chosen x = k + / R + . Eval- 
uation of the integral in equation (10) is standard and 
yields 



N 



LO 



dxdh 



2(27r) 3 x(l -x) 



^*(ar > k x )^(a :> k J -) = l ) (12) 



a simple overlap of the two-body wave function. 

To analyze the normalization to first order in G — G, 
we expand equation (ph to first order 



N 



LO 



SN ■ 



= 1, 



(13) 



where N LO is the integral appearing in ( |l2| ) and the first- 
order correction arising from (ra) is 



SN = 



2R+ 



ill 



G 



d 



dR- 



G' 



G-GWG 



GV G-G 



)( 



dR 



-G- X )G 
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The presence of G merely subtracts the leading-order re- 
sult N ho . Considering for the moment just the first term 
in the above equation (and omitting the subtraction G), 
we have 



d 4 k d 4 p 



(2tt) 4 (2tt) 4 



(l-x)7*(a;,k- L |M 2 )d(fc) 



x d(R-k) 2 V(k,p)d{p)d(R-p)j(y,p x \M 2 ). (15) 

The minus-momentum integrals above are similar to 
those considered in derivingthe bound-state equation to 
leading order in Appendix EL The only difference is the 
double pole due to the extra propagator d(R — k). With 



x = k + / R + and y = p + / R + , for x > y we avoid pick- 
ing up the residue at the double pole and the result is 
the same as in Eq. @ after using the bound-state equa- 
tion. This term is then subtracted by the G term in 
equation (|l4|). On the other hand, when x < y we pick 
up the residue at the double pole. Part of the residue 
is subtracted by the G term; the other half depends on 
dV(k,p)/dk~ . The second term in Eq. ( |l4|) is evaluated 
identically up to {k <-> p}. Now combining the two terms 
and their relevant 9 functions, we can rewrite the result 
using th e ex plicit form of the one-boson exchange poten- 
tial Eq. A3, namely 



SN 



dxdk. 



dydp A 



2(2tt) 3 x(1 - x) 2(2vr) 3 y(l - y) 



V>*(x,k J 



dM 2 



V(x,k ± ;y,p ± \M 2 ) Wp x ). 



(16) 



r 



Thus although the covariant derivative's action on the 
potential vanishes, we can manipulate the correction to 
the normalization into the form of a derivative's action on 
the light-front, time-ordered potential. With this form, 
we can compare to the familiar nonvalence probability 
discussed in Appendix [b] (in the frame where R + = R~~ = 
M/v2 with M/y/2 as the eigenvalue of the light-front 
Hamiltonian, denoted p~ in Eq. (B8)). 

Here we have seen that the normalization of the light- 
cone wave function includes effects from higher Fock 
states. Since this normalization condition stemmed from 
a diagonal matrix element of a pseudo current, we should 
not be surprised that matrix elements of the electromag- 
netic current, when treated in this reduction scheme, pick 
up contributions from higher Fock states. These higher 
Fock contributions appear explicitly and are the subject 
of the next section. 



III. APPLICATION TO GPDS 

Having worked through matrix elements of the electro- 
magnetic current in a frame where A + ^ we can 
now make the connection to generalized parton distribu- 
tions (GPDs). These distributions, which in some sense 
are the natural interpolating functions between form fac- 
tors and quark distribution functions, turn up in a variety 
of hard exclusive processes, e.g. deeply virtual Comp- 
ton scattering, wide-angle Compton scattering and the 
electro-production of mesons |l3J. The scattering ampli- 
tude for these processes factorizes into a convolution of a 
hard part (calculable from perturbative QCD) and a soft 
part which the GPDs encode. Since light-cone correla- 
tions are probed in these hard processes, the soft physics 
has a simple interpretation and expression in terms of 



In this section, we cast 
in the language of GPDs 



light-front wave functions 114 
our results for form factors J? 
and the light cone Fock space expansion. The (3 + 1)- 
dimensional expressions for form factors are presented 
in Appendix |A|. Additionally one can obtain these re- 
sults directly from time-ordered perturbation theory us- 
ing two-body projection operators as explicated in Ap- 
pendix [b| 

The GPD for our meson model is defined by a non- 
diagonal matrix element of bilocal field operators 



■in 



FU-cM = / ^e" p+ M* P ,| q(y-)id + q(0) |* P >, 

J17) 

where q(x) denotes the quark field operator and = 

_ Qt* Comparing to the current matrix element J M in 
Appendix [A|, the definition of the GPD leads immediately 
to the sum rule 



dx 



l-C/2 



F(x,C,t)=F(t). 



(18) 



Hence one can calculate these distributions from the in- 
tegrand of the form factor. In this way, the light-cone 



correlation defined in Eq. (17) has a natural description 
in terms of light-front time-ordered perturbation theory, 
e.g. for x > £ the relevant graphs contributing to the 
GPD are in Figure |^ and ||, and those for x < £ are in 
Figure |(| 



A. Continuity 

Conversion of the contributions to the form factor into 
GPDs is straightforward using Eq. (|l8|) ; we merely re- 
move —2iP + and J dx from Eqs. ( AlCj - |A15| ). In order 
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for the deeply virtual Compton scattering amplitude to 
factorize into hard and soft pieces (at leading twist), the 
GPDs F(x,(,t) must be continuous at x = £■ Main- 
taining continuity at the crossover is more pressing be- 
cause experiments which measure the beam-spin asym- 
metry are limited to the crossover [fl5[ . The leading-order 
expressions are continuous. This is easy to see since the 
contribution for x < £ is identically zero. The valence 
contribution for i > ( is a convolution of wave func- 
tions one of which is ip*(x\ . . .) which is probed at the 
end point since x' = j^- — * 0. From the bound-state 
equation Eq. (A.5), we see the two-body wave function 
vanishes quadratically at the end points. Taking into ac- 
count the overall weight x'^ 1 , the valence piece vanishes 
linearly at the crossover. At leading order then, conti- 
nuity is maintained at the crossover, while the derivative 
is discontinuous. Working only in the valence sector, va- 
lence quark models will never be of any use to beam-spin 
asymmetry measurements since the value at the crossover 
requires one wave function to be at an end point. In the 
three-body bound state problem (e.g. the nucleon), the 
valence GPD will vanish only if the three-body interac- 
tion is non-singular at the end points (which is physically 



reasonable and perturbatively true). 

Let us now check the next-to-leading order contribu- 
tions to the GPD for continuity. First we shall deal with 
the term stemming from iter ating the Bethe-Salpeter 
equation for the initial state (Ala) (see diagram B of 
Figure |^). Since there is no Z- graph generated from iter- 
ating the initial state, we expect this contribution to van- 
ish. Looking at the expression, we see again ip*(x' , . . .)/x' 
which vanishes linearly asi->(. Moreover there are the 
interaction terms: D(y, p- 1 ; x, |M 2 ) which is finite as 
x — > C, and 



D(y',p^;x',k^|A/V= 



(19) 



which vanishes at the crossover. Thus not only does the 
initial-state iteration term vanish at the crossover, its 
derivative does so as well. 



Now we investiga te th e Born terms All (see diagram 
A of Figure ||) and A12 (see diagram D of Fig ure [|) at 
the crossover. Approaching £ from above (All), we have 
the finite contribution at the crossover 



-ftp (y , p M ) ^ | ; - i N2 ; _ 2 D(y, p^; C, ^ \M )ij}{y, p ^ 



(167r3)2y(l -y)y> 



(k^+A-L) 2 



(20) 



On the other hand, approaching the crossover from below ( [A 12] ) we have to deal with singularities as x" = x/C, — ► 1. 
Writing out the propagator for the quark-antiquark pair heading off to annihilation, we see 



DwCAk'^l*)^- 



1 - x" 



(k 1 - +A^) 2 + m 2 



(21) 



This linear vanishing cancels the weight (1 — x")^ 1 . Taking the limit x — > C then produces equation (|20| ) and thus 
the Born terms are continuous. 

Lastly we must see how the final-state iteration terms match up at the crossover. Using Eq. ( A14 ) to approach £ 
from above (see diagram C of Figure ^|), we have the contribution 



final 



(22) 



Approac hing £ from below (see diagram E of Figure ^J), we utilize equation ( pl| ) in taking the limit of (A15). The 
result is ( |22|) and hence we have demonstrated continuity to first order, i.e. 



F(Ct ) t)=F((,(;,t) B °™ + F(( ) (,t) 



final 



(23) 



no matter how we approach x = (. 



B. Fock space representation 

We now write the GPDs in terms of Fock component overlaps. In the diagonal overlap region x > £ this will be a 
mere rewriting of our results, while there is a subtlety for the non-diagonal overlaps. To handle the zeroth-order term 
( |A1C| ), we define the two-body Fock component as 

MxiM^zM) = -^=il>(xx,)s^), (24) 

y/XlX2 
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noting that the relative transverse momentum can be defined as k^j = rz^kj 1 — xik^-. In terms of Eq. 
appears as 



LO 



where the primed variables are given by 

J x i = T=*c 

1 L-Zi. 



J'=l,2 



the GPD 



(25) 



j — i-c 

-L - 1,1- 



k+ - ^A X , for t ^ j [k^ = k+ + (1 - x' J )A J 
and the integration measure is given by 



(26) 



N 



N 



[cIx]n = Y\dxj S[ l - 



i=l 



W TV 



(27) 
(28) 



Notice the sum over transverse momenta in the delta function is zero since our initial meson has P 1 " = 0. The sum 
over j in Eq. ( ^5| ) produces the overall factor of two for our case of equally massive (equally charged) constituents. 

To cast the next-to-leading order expressions for x > £ in terms of diagonal Fock space overlaps, we must write 
out the three-body Fock component. Looking at the diagrams in Figure ^[ it is constructed from the two-body wave 
function 



tp 3 (xi,k^) = g 



2(2tt) 3 

y/XiX 2 X 3 



[dyh{dp J 



0(yi - x 1 )x 3 8{y 2 - z 3 )<KP2 ~ H)D(yi,pi;xiM\M 2 ) 



+ %2 - x a )x 1 S(yi - xt)S(pi - ki)D{y 2l p^; x 3 , k^|M 2 ) 



, (29) 



where i runs from one to three and the label j, which stems from the integration measure, runs from one to two. We 
discuss how to obtain this three-body wave function directly from time-ordered perturbation theory in Appendix |^. 
Using ip 3 in Eq. (E9), the terms in the GPD at first order in the weak coupling can then be written compactly as 



F(x,C,t) 



NLO 



'(X ~ C) 

i-C 



idxUdk^ J2 6(x- Xj )r 3 (x^) 

.3=1,3 



XjXj 



i/) 3 (xi,k±). 



(30) 



One can verify that the diagrams in Figure |^ are gener- 
ated by (|3C)|). Additionally there is a fourth diagram gen- 
erated by Eq. ( pp[ ) which does not appear in the figure. 
This missing diagram is characterized by the spectator 
quark's one-loop self interaction and is absent since we 
have ignored f(k 2 ) and the scale dependence of light-cone 
wave functions. 2 The absence of this diagram does not 
affect continuity at the crossover. The missing diagram 
vanishes at x = £ since the final-state wave function is 
ip*(x', . . .). Lastly we note the above Fock component 
overlaps satisfy the positivity constraint for a composite 



2 Recently the scale dependence of light-cone Fock components has 
been investigated Jig] . 



scalar composed of scalar constituents |L7j ■ 

Now we must come to terms with the non-diagonal 
overlap region, x < (. At first order, the diagrams of Fig- 
ure^ correspond to four-to-two Fock component overlaps. 
We have been cava lier a bout time ordering, however. The 
expressions Eqs. A12 and A15 do not correspond to 



time-ordered graphs. Both terms contain a product of 
time-ordered propagators: one for the two quarks lead- 
ing to the final-state vertex and another for the quark- 
antiquark pair heading off to annihilation. But for an 
interpretation in terms of a four-body wave function, all 
four particles must propagate at the same time. This is 
a subtle issue as a graph containing the product of two 
independently time-ordered pieces (where one leads to 
a bound-state vertex) corresponds to a sum of infinitely 
many time-ordered graphs. It is easiest to write out the 
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terms of concern in terms of the propagators' poles. The 
quark, anti-quark heading to annihilation have propaga- 
tors d(k) and d{k + A) and poles we label k~ and k~, 
respectively. The remaining propagators of interest d(p') 
and d(P' — p') have poles p~ and , respectively. We 
can then manipulate as follows 



1 



1 



l ( l l 



Pb P a ^ c fca \ P Pa k c k a 



(31) 



In this form, we have produced the correct energy de- 
nominator for the instant of light-front time where four 
particles are propagating. Multiplying this denomina- 
tor by the three-body wave function yields the four-body 
wave function (up to constants). This is the part of the 
four-body wave function relevant for GPDs (there are ad- 
ditional pieces for two-quark, two-boson states, see Ap- 
pendix |b|). In the resulting sum (31), the first term will 
produce the two-body wave function for the final state 
and we will have a genuine four-to-two overlap. We do 
not write this out explicitly. 

The second term in Eq. (|3l|), however, contains again 
the propagator for the pair heading to annihilation. Us- 
ing the light-front Bethe-Salpeter equation for the vertex 
(which contains infinitely many times) we can introduce 
a factor of the time-ordered interaction. The resulting 
product of independent time ordcrings can again be ma- 
nipulated as in Eq. (|3l|). The result produces another 
overall four-body denominator which contributes to the 
four-body Fock component of the initial state. Since we 
iterated the interaction, however, this new contribution 
is no longer at leading order and can be neglected. Thus 
the second term in (|3l]) does not contribute at this order. 

Having manipulated the GPDs into non-diagonal over- 
laps for x < C, we must wonder if continuity at the 
crossover is still maintained. In the limit x — > £ the 
light-front energy of the struck quark goes to infinity. 
Consequently k~ , which contains this on shell energy, is 
infinite and dominates the four-body energy denomina- 
tor. This is identical to the reasoning in Eqs. A12| and 
A15 where instead of the four-body denominator, we have 
D\v(x", k +x"A \t) which is dominated by k~ at the 
crossover. Either way, we arrive at the expressions found 
above for the crossover ( pp| ) and p2]). 

Having cast our expressions for generalized parton dis- 
tributions in terms of the Fock components, we can en- 
large our understanding of the sum rule and continuity 
at the crossover. Both must deal with the relation be- 
tween higher Fock components. The way (^-dependence 
disappears from (|l^ ) mandates a relation between the di- 
agonal and non-diagonal Fock component overlaps that 
make up the GPD. The relation between Fock compo- 
nents must follow from the field-theoretic equations of 
motion. Continuity itself is a special case of the rela- 
tion between Fock components, specifically at the end 




FIG. 1: The triangle diagram for the time-like, pion form 
factor 



points. Above we have seen our expressions are contin- 
uous (and non-vanishing) at the crossover and explicitly 
that the three- and four-body components match at the 
end point (where x— £ = 0). This weak binding model for 
behavior at the crossover is a simple example of the rela- 
tions between Fock components at the end points (see 
also p8| ). More general relations must be permitted 
from the equations of motion to guarantee Lorentz co- 
variance (e.g. in the structure of the Mellin moments 
of the GPDs, of which the sum rule is a special case). 
Here, of course, Lorentz symmetry is broken. Infinitely 
many light-cone time-ordered graphs are needed in the 
reduced kernel to reproduce the covariant one-boson ex- 
change (|^). Thus exactly satisfying polynomiality re- 
quires not only infinitely many exchanges in the kernel 
but contributions from infinitely many Fock components. 
It should be possible, however, to show how the sum rule 
and polynomiality are improved order by order. 



IV. APPLICATION TO GDAS 

Below we study time-like form factors to demonstrate 
the versatility of this approach and make the connec- 
tion to the generalized distribution amplitudes (GDAs) 
for this model. Analogous to GPDs, GDAs encode the 
soft physics of two-meson production and can thus be 
thought of as crossed versions of the GPDs. The GDAs 
enter in convolutions for various two-meson production 
amplitudes ]I9|] . These distribution functions as well as 
time-like form factors are a theoretical challenge for light- 
front dynamics, since there is no direct decomposition in 
terms of meson Fock components alone. Furthermore, 
we shall see the leading-order expressions are non- valence 
contributions (which necessarily excludes a description in 
terms of most constituent quark models). 

The time-like form factor F(s) for our model meson is 
defined by (see Figure |1) 



(% *p/|r" \o) = -i(p-j/)»F( 8 ) 



(32) 
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FIG. 2: Leading order diagrams for the time-like, pion form factor. 



where s = (p + p') 2 is the center of mass energy squared. 
Now define P» = p^ + p' and C = P + / p+ - We can 
work out the kinematics of this reaction in a frame where 
P ± = 

P- = 



V 



2P+ 

(1-0* 
2P+ 

*(i-0C 



M 



where M is the meson mass. 

Similar to GPDs, the GDA for our model has a defini- 
tion in terms of a non-diagonal matrix element of bilocal 
field operators 



<r>u.,.s,= / ^e- p+ - 



(* p * p ,\q(x-)i d + q(Q) |0). 

(34) 

Such a definition of the GDA leads directly to a sum rule 
for the time-like form factor- 



where we have made use of z = k + / P + , z' = and 
z" = z/(. Recall ^(x\M 2 ) cx 9[x(l - x)]. This translates 
to: < z" < 1 and < z' < 1, and hence we do not pick 
up a contribution at zeroth order in the coupling. 



To work at first order, we pick up three terms anal- 
ogous to those in Appendix [A|. We denote these as 
5Jif,SJ£ and Sj£,. The Born term for the three-point 
(33) electromagnetic vertex SJ^ is quite simple. For the 
same reason as the zeroth-order result, the restriction 
of r ){x\M 2 ) oc 9[x(l — x)] and momentum conservation 
force the contribution 8Jti to vanish. This leaves us to 
consider only diagrams that arise from iteration of the 
Bcthc-Salpctcr equation of cither final-state meson. 



Considering first the term 5J£, we have the contribu- 
tion to the GDA 



dz 



2C-1 



$(z,C,s) = F(s), 



(35) 



and hence a means to calculate $ from the integrand of 
the time-like form factor. 

Taking the appropriate residues of the five-point func- 
tion, we arrive at Fig. |l| for the time-like form factor. 
Keeping only the leading order piece of the electromag- 
netic vertex r M , we have 



$(z,C,s) = iP+(2z-l) 



dk-dk 1 



x 7 *(z",k ± -z'V \M 2 )G(k,p)d- 1 (k -p) 

x G{P - k, p')7* (z, -(1 - z') p- 1 |M 2 ), (36) 



$ p (z,C,s) = zP+(2z-l) 



dk-dk 1 - d 4 q 



(2tt) 4 (2tt) 4 

x 7* (y", q X ~ y" p X )G(q,p)V(q, k)G(k,p)G(P - k, p') 

x d-\k-p)r(z'M -(1 - z') p-L), (37) 

where we have chosen to abbreviate y = q + / P + and 
hence the label y" = y/Q. We have customarily omitted 
the subtracted term containing G, which is zero because 
there is no leading-order term to subtract. Requiring 
wave function vertices mandates < y" < 1 and C, < z < 
1. 

Thus $ p produces one contribution to the GDA 



(I6n d yc J z{l-z) z'y"(l - y") 



x D(z, k x ; y, q x | S )V*(z', k x -(1 - z') p^* (y", q X - y" p X ). (38) 
As a contribution to the time-like form factor, we can interpret Eq. ( |3"§| ) as the time-ordered diagram 6 of Figure ^| 
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At first order in the weak coupling, we have one term remaining to consider 5 j£, . Again omitting the superfluous 
subtraction of G, we have 

dV(^C,s)^P + (2z-l) |^^^7*(Ak- L -2 ,/ p- L |Af a )G(*, J ,)d- 1 (p-*) 

x G(P - k, p')V(P -k,P- q)G(P - q, P ')7*(y', -(1 - /) P X l^ 2 ), (39) 
with y' = jE§ • Both y' and z are restricted: < y' < 1 and < z < £, and the final contribution to the GDA is 

(l&iT A yC J (1 - z)z"(l - z")y'(l - y'J y-z 

x D(y, q^; 2, | S )^(z", k^ - z" p^V, ^ + y' P^), (40) 



In this form we recognize this contribution as diagram a 
of Figure ||. 

Having found the leading non-vanishing contribution 
to the GDA namely $ = <£> p + <£y , we observe that the 



higher Fock components derived in section III (as well 
as in Appendix |b|) do not fit naturally into (3^) or (40). 
One needs a Fock space expansion for the photon wave 
function in order to have an expression for the GDA in 
terms of various Fock component overlaps. With the ex- 
pressions derived for the GDA we can use Eq. (|35|) to 
obtain the time-like form factor. 



V. SUMMARY 

Above we have investigated various applications of the 
light- front reduction of current matrix elements. First we 
considered the normalization of the light-front wave func- 
tion in the reduction formalism deriving Eq. (||) from the 
covariant normalization Eq. (0) . The complicated form of 
the reduced normalization was linked to effects of higher 
Fock components (which we illustrated by using the lad- 
der model (^|) in perturbation theory). Using the explicit 
form of the leading-order kernel, we were able to derive 
Eq. ([l6]) , which is the familiar many-body normalization 
condition (|B8|). 

In Appendix |A|, we reviewed the derivation of the form 
factor at next-to- leading order in the (3 + l)-dimensional 
ladder model. These expressions were then converted 
into the GPD for the model. Continuity of these dis- 
tributions at the crossover (where the plus momentum 
of the struck quark is equal to the plus component of 
the momentum transfer) was explicitly demonstrated, cf 
Eq. ((2^). Connection was made to the overlap repre- 
sentation of GPDs by constructing the three-body wave 
function to leading order in perturbation theory. As 
a check on our results, we also reviewed the construc- 
tion of higher Fock states from the valence sector in old- 
fashioned time-ordered perturbation theory (Appendix 
|b|). The derived overlaps Eqs. ( p5| ) and ( |30| ) satisfy the 
relevant positivity constraint. The non- vanishing of the 



GPDs at the crossover could then be tied to higher Fock 
components, specifically at vanishing plus momentum, 
and are hence essential for any phenomenological mod- 
eling of these distributions. This rewriting allowed us 
to understand how continuity arises perturbatively from 
the small- a; behavior of Fock state wave functions. In 
perturbation theory, the diagonal valence overlap van- 
ishes at the crossover, while the higher Fock component 
overlaps do not. In general the n-to-n overlap matches 
up with the (n + l)-to-(n — 1) overlap at the crossover 
due to the dominance of the rebounding quark's infinite 
energy. The same seems to be true perturbatively for a 
three-body bound state due to the nature of the kernel. 

Unfortunately issues involving Lorentz invariance 
(such as the sum rule for the electromagnetic form factor 
and the polynomiality constraints) are left untouched. 
To maintain covariance one would need infinitely many 
time-ordered exchanges in the kernel as well as infinitely 
many Fock components. It should be possible, however, 
to understand perturbatively how the ^-dependence dis- 
appears from Eq. (|l8|). This requires further relations 
between Fock components and these should be afforded 
by the field-theoretic equations of motion. 

Lastly we considered application of the reduction for- 
malism for currents to the time-like form factor. We did 
so by calculating the ladder model's GDA Eqs. (Hf-[Io|), 
which is related to the time-like form factor via the sum 
rule in Eq. (|3^), systematically in perturbation theory. 
This is in contrast to the non-existent Fock space expan- 
sion for these types of processes. 

With the formalism explored here, one could use phe- 
nomenological Lagrangian based models to explore both 
generalized parton distributions and generalized distribu- 
tions amplitudes within the light-front framework. Such 
an investigation is interesting not only for testing phe- 
nomenological models, but also for anticipating problems 
for approximate non-perturbative solutions for the light- 
cone Fock states. Nonetheless more model studies are 
warranted before truly realistic calculations can be pur- 
sued. 
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+ 



FIG. 3: Diagrammatic representation of the one-boson ex- 
change potential V appearing in Eq. (A3). 
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APPENDIX A: WAVE FUNCTIONS AND FORM 
FACTORS IN (3 + 1) DIMENSIONS 

In this Appendix we collect results relevant above for 
wave functions and form factors in the (3+l)-dimensional 
ladder model Eq. <M). 




k + A 



FIG. 4: The leading-order diagram for the electromagnetic 
form factor 



For reference, the bound-state equation (Al) appears 



as 



y) 



xV{x,\^;y^\M')i,{ y , V ^) (A5) 



1. Wave functions 

Using the Bethe-Salpeter equation (||) and the defini- 
tion of the light-cone wave function \4>r) @) we have the 
light-cone bound-state equation 

= g(R)w(R)\ip R ), (Al) 

where w(R) is the reduced auxiliary kernel 

w(R) = g-\R) G(R)W(R)G(R) g^(R), (A2) 

with W(R) defined in Eq. (§). 

To leading order in G — G, one calculates w(R) for the 
ladder model to be: 



Form factors 



V(x,k^:y,p^\M z ) 

g 2 

x-y 



2\ — 



xR + ,^ | w(R) | yi? + ,p ± ) 
(x-y)D(x,k ± ;y,P ± \M 2 ) 
- {(x, k x ) <— (y, p^)}] 9[x(l - x)]6[y(l - y)}, (A3) 
where we have defined 



D-\x,k x ;y,p x \M 2 )=M 2 



(k^-p^) 2 + M 2 

x-y 



1-x 



(A4) 



and taken R- 1 = 0. Graphically this one-boson exchange 
potential Eq. (A3) is depicted in Figure [| 



To calculate form factors, we use the electromagnetic 
vertex constructed in [Q up to the first Born approx- 
imation (notice that the ladder model's gauged interac- 
tion V = 0) 



(5 



(A6) 



where 9 M denotes the electromagnetic coupling to 
scalars. 

Now using Eq. @ to first order in G — G, the matrix 
element = {*S> p^T^^-A)^ P ) then appears 



J"« ( 7P ,| G(P'){l + V(P')(G(P')-G(P')) 
x (d^{-A)d^ + V(-A)G(-A)d^(-A)d^ 
x (i + (G(P)-G(P))V(P))g(P)\ | 7p ) 
= (j£o + *J? + + K) + o[v% 
with the leading-order result 

■tfo = (tH \G(P')d^(-A)d^G{P)\ | 7 p). 
The first-order terms are 

SJt = M G(P')^(-A)rf 2 - 1 



(A7) 



(A8) 



(g(P)-G(P))v(P)G(P)\ | 7 p) 
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FIG. 5: Diagrams which contribute to the form factor at next-to-leading order (for x > £). 




FIG. 6: The remaining diagrams (characterized by x < Q for the electromagnetic form factor at next-to-leading order. 



6J» = { 7P ,\ G(P')V(P') 



x (G(P')-G(P'))d^-A)d^G(P)\ | 7P ) 

sji; = { lP ,\ |g(p')(f(-a)g(-a)^(-a)) 

xd^G(P) \<y P ). (A9) 



As outlined in ffi], Eqs. ( A.8) and (A.S) can be evaluated 
by residues being careful to remove two-particle reducible 
contributions by utilizing ( |A5| ) . Here we state the results 
of these calculations in (3 + 1) dimensions. We denote 
A M as the momentum transfer and define A + = — (P + 
(see Figure 0). The leading-order result appears 



Jlo = -2^ 4 



2{2tt) 6 x{1 — x)x' 



where x' = jz^ and k /J - = k^ +(1 — x 1 ) A 1 - denotes the momentum of the final state. Using J M 
Eq. ([AH]) reduces to the Droll- Yan formula |lj for C = 0. 

The first of the leading order corrections is the Born term <5J+. For a; > ( we have 



(A10) 

-i(P + P r YF(t), 



+ _ +2iP+ r 9(x-C)dxdk ± dydp ± (2x-C) 
7 (»>0 ~ (16tt 3 ) 2 J xx'y{\ - y)y' 



x r (y\ p' x )D(y', p' x ; x', k' x ]M 2 ) gH{v _ x) D(y, p^; x, \M 2 )^(y, p x ), (All) 

y 

where y' — jz^ and p /J - = p- 1 +(1 — y') A- 1 . This contribution corresponds to diagram A in Figure On the other 
hand, for a; < £ we have 

+ +2iP+ r 9(C-x)dxdk ± dydp^(2x-C)/C 

7 (*<C) ~ (16tt 3 ) 2 7 y(l-y)y'x"(l-x") 

x ^(y',p' ± )D w (x",k" ± |t) g2% ~ ^ Dd/.p^ig.k 1 \M 2 )i,{y,p^), (A12) 

where x" = x/C, and k" 1 - = k- 1 +x" A x denotes the photon's relative momentum. This expression is diagram D in 
Figure @. 
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The next leading-order term is the initial-state iteration 5J^~ . The only contribution is for x > £, namely 
+2i P+ f 6{x- Qdxd k- 1 dyd p- 1 (2x - C) 



SX 



(16tt 3 ) 2 



xx'(l — x')y(l — y) 



x^(^,k^) J D(y',p' J -;x',k' J -|M^) 



Z^p^x^lM 2 )^^), (A13) 



which corresponds to diagram B of Figure |. 

Lastly there is the final-state iteration term 8Jj . For x > ( we have 



+ _ +2iP+ f e(x-C)dxdk- L dy'dp' ± (2x~Q 
f (*>0 " (16tt 3 ) 2 



x(l-x)x'y'(l-y') 



xf(y',p' 1 )B(y',p' 1 ;x',k' 1 |A/ 



^ g 2 9(y-x) 
y-x 



D{y^ ]X ^\M % )^{x,\^), (A14) 



where implicitly y = Q + (1 — £) y' and p^ = p' x — (1 — y') A^. This corresponds to diagram C in Figure |J. While 
for x < £, the expression 



_l_ +2^P+ f 6(C - xjdxdk 1 - dy' dp ,± {2x - Q/C 

/ (x<c) - (i 67r 3)2 _/ (l-x) x"(l-x")y'(l-y') 

x^(y , ,P ,± )Dw(x",k" i |i) 



gH^y-x) 
y-x 



D^p^^lM 2 )^,^), (A15) 



corresponds to diagram _E of Figure |[ Eqs. (A1C-A15) are then the complete expressions for the form factor up to 
first order. 



APPENDIX B: OLD-FASHIONED 
TIME-ORDERED PERTURBATION THEORY 

The results of this paper can similarly be achieved 
directly from "old-fashioned" time-ordered perturbation 
theory in a form which utilizes projecting onto the two- 
body subspace of the full Fock space. For a nice, com- 
plete discussion of this formalism for the light-cone ladder 
model, see ||(J. In this Appendix, we show how to derive 
higher Fock space components in this formalism, thereby 
demonstrating the generation of higher components from 
the lowest sector we found indirectly for GPDs and form 



factors in section III 



We write the light-cone Hamiltonian as a sum of a free 
piece and an interacting piece which carries an explicit 
power of the weak coupling g. In an obvious notation 
this is 



P- = P~ 



9P7 



(Bl) 



The free term P~ is diagonal in the Fock state basis, 
while the interaction generally mixes components of dif- 
ferent particle number (in the scalar model we consider 
above, the interaction is completely off-diagonal since 
there are no instantaneous terms). Let us suppose that 
in the full Fock basis, we have an eigenstate of the Hamil- 
tonian, i.e. 



(B2) 



P0+9P7) W) = v 



where the eigenvalue is labeled by p~. Since the coupling 
is presumed small, the mixing of Fock components with 
a large number of particles will be small. Thus one imag- 
ines our bound state will be dominated by the two-body 
Fock component. 

To make this observation formal, we define projections 
operators on the Fock space V and Q in the usual sense. 
The operator V projects out only the two-particle sub- 
space of the full Fock space and hence Q projects out the 
compliment. Let us define the action of these operators 
on our eigenstate 



(B3) 
(B4) 



As is well known, combination of Eqs. (B3) and (B4) 
leads to the following equation for the two-body Fock 
component 



i^hfe> = P"hfe>, 

where the effective two-body Hamiltonian is 



P 



vv 



^eff 



= P. 



VV 



p 



1 



V Q , 



-P, 



QQ 



QV> 



(B5) 

(B6) 
(B7) 



and we have defined the following notation for any opera- 
tors A and B, P^ B = AP~B. The effective two-body in- 
teraction V c f{ defined in equation ( B6 ) is dependent upon 
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the energy eigenvalue p~ since we have suppressed the de- 
grees of freedom of the Q subspace. The relation between 
the Q-space probability (i.e. the non-valence contribu- 
tion) and the effective interaction appears 



(V'Q 



d 



dp 



- (^2 | Ves\lp2 



(B8) 



In a weak-binding limit, we can series expand the ef- 
fective interaction in powers of the coupling and thereby 
re-derive the light-front potential. Given that every bo- 
son emitted must be absorbed in the two-quark sector, 
we can have only an even number of interactions and 
hence 



V cS = VgPf 




(B9) 

So, for example, at leading order we have all possible 
ways to propagate from the two-body sector and back 
with only two interactions in between. The diagrams in 
Figure || correspond to the two possibilities distinguished 

by the action of — between interactions. At the next 

p~ —Po 

order, we have all possible ways to propagate from two 
bodies to two bodies with four interactions in between, 
etc. 

To generate higher Fock components from the two- 
body sector, we necessaril y m ust l ook at the Q-space 
state arrived at from Eqs. (B3) and (B4) 



I^Q> = 



1 



p- - QV 
Q Q 



(BIO) 



To generate an n-body Fock component from this state, 
we merely act with an n-body projection operator which 
we shall denote Q n . Similar to the above, we expand in 
powers of the coupling to find 



\lpn) = Qn 



p- -Po 



n=0 



-Po 



QgPf IV> 2 > 



(Bll) 

For example, the leading-order three-body state is ob- 
tained by attaching a boson to a quark line in the only 
two possible ways (and adding the light-front energy de- 
nominator at the end). With these three-body states, 
we can consider all possible three-to-three overlaps that 
would contribute to the form factor. These are de- 
picted in Figure [| (with the exception of a quark self- 
interaction). The four-body sector is richer since there 
are two-boson, two-quark states as well as four-quark 
states. The two-to-four overlaps required for GPDs must 
have four quarks. At leading order, we generate the 
diagrams encountered above in Figure | Not surpris- 
ingly directly applying time-ordered perturbation theory 
from a light-front Hamiltonian agrees with our derivation 
above from covariant perturbation theory in the Bcthc- 
Salpeter formalism. 
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